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Psj . Abstract: In order to deal with a large cosmological constant a relaxation mechanism 

, based on modified gravity has been proposed recently. By virtue of this mechanism the 

ff^ \ effect of the vacuum energy density of a given quantum field/string theory (no matter how 

big is its initial value in the early universe) can be neutralized dynamically, i.e. without fine 
tuning, and hence a Big Bang-like evolution of the cosmos becomes possible. Remarkably, 



0\ 

■ a large class {-^^1 of models of this kind, namely capable of dynamically adjusting the 

vacuum energy irrespective of its value and size, has been identified. In this paper, we 
carefully put them to the experimental test. By performing a joint likelihood analysis 
we confront these models with the most recent observational data on type la supernovae 
(SNIa), the Cosmic Microwave Background (CMB), the Baryonic Acoustic Oscillations 
\ (BAO) and the high redshift data on the expansion rate, so as to determine which ones 

are the most favored by observations. We compare the optimal relaxation models 
found by this method with the standard or concordance ACDM model, and find that some 
of these models may appear as almost indistinguishable from it. Interestingly enough, 
this shows that it is possible to construct viable solutions to the tough cosmological fine 
tuning problem with models that display the same basic phenomenological features as the 
concordance model. 
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1. Introduction 

Current observations [0, |2| indicate that the cosmological constant (CC) A, or equivalently 
the vacuum energy density p\ = h/ (SttGn), is non- vanishing and positive, and of the order 
of ~ 10"^'' GeV^. This value is close to 70% of the present critical energy density in 
our universe, p^, and therefore is very small in Particle Physics units (being equivalent to 
having a mass density of a few protons per cubic meter). In itself, the inclusion of the tiny 
value p^ as a cosmological term in the gravity action would not be a problem. However, 
the modern fundamental physical theories suggest that large contributions to the classical 
CC parameter (in particular those sourced by quantum fluctuations of the matter fields) 
are induced, resulting into a huge initial CC value of order ~ Mj^ which is fed into the 
cosmos already in the phase transitions of the early universe, with Mx in the range between 
the electro-weak (EW) scale ~ 100 GeV and the Planck scale ~ 10^*^ GeV. Even if "only" 
the EW scale of the Standard Model (SM) of Particle Physics would be involved in the 
induced CC value, the discrepancy with respect to the observed value p\ entails 55 orders 
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of magnitude^. This is of course preposterous. Despite none of these vacuum contributions 
is directly supported by experimental evidence yet, there is hardly a theory beyond the 
SM that does not come with a large value of A = Stt Gn Pa induced by the huge vacuum 
energy density p\ predicted by the theory and whose origin may be both from the zero point 
vacuum fluctuations of the matter fields as well as from the spontaneous symmetry breaking 
of the gauge theories. In principle, all models with scalar fields (in particular the SM with 
its Higgs boson sector) end up with a large vacuum energy density p\ associated to their 
potentials. For example, from the SM Higgs potential we expect p\ ~ 10^ GeV^, which 
explains the aforementioned 55 orders of magnitude discrepancy and corresponding 55*^ 
digit fine tuning CC problem. The huge hierarchy between the predicted and the observed 
CC, the so-called "old CC problem" Q, is one of the biggest mysteries of theoretical physics 
of all times. To avoid removing the huge initial value of the vacuum energy density p\ 
by hand with an extremely fine tuned counterterm ^ (giving the theory a very unpleasant 
appearance), alternative dark energy models have been suggested. Typically, these models 
replace the large CC by a dynamical energy source [see [|[ and references therein]. This 
seems reasonable for describing the observed late-time acceleration, but generally it still 
requires that the large CC has been removed somehow. In the end, it means that a fine 
tuned counterterm has been tacitly assumed |^] . 

It is encouraging that time-varying vacuum energy models inspired by the principles 
of QFT can be suggested and may hopefully provide an alternative and more efficient 
explanation for the dynamical nature of the vacuum energy . Recently the analysis of 
some well motivated models of this kind versus the observations has shown that their 
phenomenological status is perfectly reasonable and promising [§, |9| - see also the recent 
works |1C1| ]. However, in order to solve the old CC problem in this context, we need to 
focus on models that dynamically counteract the large initial p\. The idea, in a nutshell, 
is: 1) to obtain an effective PKes{t) (the measured one) which satisfies |/OAefr(OI \pW 
at all times t and without fine-tuning; 2) to insure that PKes{t) preserves the standard 
cosmic evolution (i.e. the correct sequence of radiation, matter and dark energy dominated 
epochs); and finally 3) to make sure that PKes{t) is able to reproduce the measured value 
of the vacuum energy at present, /5Aeff(io) ~ p\-, again without fine-tuning. Only in this 
way we can guarantee a low-curvature universe similar to ours, namely R = 8tt G^PAeS ~ 
SttGnPc ~ Hq, with Ho ~ 10~^^ GeV the current Hubble rate. The question of course is: 
is that program really possible? 

Obviously, in order to follow this road in a successful way it is necessary to go beyond 
the common dark energy models and face a class of scenarios where the large CC is not 
hidden under the rug, so to speak, but rather it is permanently considered as a fundamental 
ingredient of the overall theory of the cosmos. In this work, we analyze a class (probably not 
unique) of models operating along this line, namely the CC relaxation mechanism based on 

^Let us note that the QCD scale of the strong interactions, Aqcd ~ 100 MeV, could also trigger an 
induced vacuum energy density pA ~ Aqcd) which would then be roughly 43 orders of magnitude larger 
than We have nevertheless taken the EW value because it is larger and is considered also as a robust 
fundamental scale in the structure of the SM of Particle Physics. 

^See e.g. a detailed account of the old fine tuning CC problem in Sect. 2, and Appendix B, of Ref. fl]. 
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F{R,Q) modified gravity 11, These models are intimately related with the family 
of the so-called AXCDM models of the cosmic evolution |jl^ (in which X is generally not 
a scalar field, but an effective quantity in the equations of motion which derives from the 
complete structure of the effective action) . They might also be connected with mechanisms 
based on matter with an inhomogeneous equation of state [0]. A first modified gravity 
approach implementing the relaxation mechanism was presented in |15]. Furthermore, a 
related model in the alternative Palatini formalism has been studied [IC], too. 

As stated, the class of models considered here is probably just a subset of a larger class 
of theories that can produce dynamical relaxation of the vacuum energy. Our work should 
ultimately be viewed as aiming at illustrative purposes only, i.e. as providing a proof of 
existence that one can construct explicitly a relaxation mechanism, if only in a moderately 
realistic form^. We cannot exclude that more sophisticated and efficient mechanisms can 
be eventually discovered which are much more realistic, but at the moment the F{R, Q) 
models serve quite well our illustrative purposes. They constitute a potentially important 
step in the long fighting of Theoretical Physics against the tough CC problem, especially 
in regard to the appalling fine tuning conundrum inherent to it. Additional work recently 
addressing the CC problem from various perspectives can be found e.g. in Refs. p8| ]. 

In contrast to late-time gravity modifications, in the relaxation mechanism an effective 
energy density prit) is permanently induced by the F{R,Q) modification to the standard 
gravity action, which counteracts dynamically (at any time of the postinflationary cosmic 
history) the effect of the large vacuum energy density p\ and provides the net value pAeff = 
P\+Pf (presumably close to the observationally measured one p^). As a result the universe 
has an expansion history similar, but not identical, to the concordance ACDM model. 
Therefore, it is important to investigate how the relaxation models perform with respect 
to observational constraints coming from the Cosmic Microwave Background (CMB) Q, 
type la supernovae (SNIa) 0], the Baryonic Acoustic Oscillations (BAO) |jl^ and the high 
redshift data on the expansion rate . 

The paper is organized as follows: in Sect. ^ we discuss different aspects of the working 
principle of the CC relaxation mechanism in modified gravity. In Sect. ^ we present the 
numerical tools that we use to perform the statistical analysis, and in Sect. Qwe apply them 
to determine the optimal relaxation model candidates with a further insight in the details 
of the relaxation mechanism. In Sect. |5| we present predictions for the redshift evolution 
of the deceleration parameter and the effective equation of state, and compare with the 
ACDM. In the last section we deliver our conclusions. Finally, in an appendix we provide 
some useful formulae discussed in the text. 



2. Vacuum relaxation in F{R, Q) modified gravity 

The general form of the complete effective action of the cosmological model in terms of the 

^Remember that in the past the dynamical adjustment of the CC was attempted through scalar field 
models jl^], but later on a general no-go theorem was formulated against them Q. Fortunately, this theorem 
can be circumvented by the relaxation mechanism, see [M for details. 
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Ricci scalar R and the Gauss-Bonnet invariant Q is given by 



/ 



1 



R- p\-F{R,Q)+C^ 



(2.1) 



where denotes the Lagrangian of the matter fields. This action represents the Einstein- 
Hilbert action extended by the term F{R, Q) defining the modification of gravity. Moreover, 
we include all vacuum energy contributions from the matter sector in the large CC term p\. 

From the above action the modified Einstein equations are derived by the variational 
principle dS/dg""^ = 0. They read 



(2.2) 



with the Einstein tensor G"^ = i?";, — cosmological term and the energy- 

momentum tensor of matter emerging from C^. 

We describe the cosmological background by the spatially flat FLRW line element 
ds^ = dt'^ — a^{t)dx'^ with the scale factor a{t) and the cosmological time t. Accordingly, 
the curvature invariants R = 6if^(l — q) and Q = —2AH'^q can be expressed in terms of the 
Hubble expansion rate H = a/a and the deceleration parameter q = —da/a^. The tensor 
components in are given by G% = -SiJ^, G'^ = -6'j{2H + SH"^) and 



1 r 
2 



F - 6{H + H'^)F^ + 6HF^ - 2AH'^{H + H'^)F^ + 2AH^F^ 
F -2{H + 3H^)F^ + AHF^ + 2F^ 
2AH'^{H + H'^)F^ + 16H{H + H^)F^ + 8H^F^ 



(2.3) 



(2.4) 



where F^ = dF/dY stand for the partial derivatives of F with respect to y = R,G- The 
F-term induces the effective energy density pp = 2E^q and pressure pp = —^E-\ implying 
that the whole effective dark energy density and pressure consist of two parts each: 



PAcS (t) = p\ + pF{t), PAefT (*) = "Pa + (^ ) • 



(2.5) 



As announced, they follow the structure of the AXCDM model |13], as both are the sum 
of a cosmological term energy density (resp. pressure) and an extra contribution whose 
structure derives from the presence of new terms in the effective action. The local energy 
density pp constitutes the "induced DE density" for this model. It is indeed induced by 



the new term of the modified gravity action (2.1). However, the only measurable DE 
density in our model is the effective vacuum energy /OAefr in (^), i.e. the sum of the initial 
(arbitrarily large) vacuum energy density p\ and the induced DE density pp- None of 
the latter, though, is individually measurable. Let us note that the induced DE density is 
covariantly self-conserved, i.e. independently of matter (which is also conserved). Indeed, 
the Bianchi identity on the FLRW background leads to the local covariant conservation 
laws 

Pn + 3H{pn+Pn)=0, (2.6) 
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which are vahd for all the individual components {n = radiation, matter and vacuum) 
with energy density pn and pressure pn- This leads immediately to the usual expressions 
Pm = P%iO'~'^ and pr = p^a~^ for cold matter and radiation respectively. As for the effective 
vacuum energy density PAefr, since p\ is a true CC term and satisfies p\ = —p\, it follows 
that the i^'-term density is self-conserved: 

PF + 3H{pF + Pf) = Pf + ^H [1+ UJF{t)] PF = 0, (2.7) 

where ujf = Wi?(t) is the corresponding (non-trivial) effective equation of state (EoS) of 
the F-term. We shall come back to it in more detail in Sect. ^ 

The cosmological evolution is therefore completely determined by the generalized Fried- 
mann equation of our model, 

3H^ ■ 

= Pni+ PAcS = Pm + Pk + PF, (2.8) 

where pm oc denotes the energy density of the (self-conserved) matter, which is mostly 
dust-like matter for the purpose of this paper (as we focus mainly on the matter dominated 
and DE epochs). 

Next, we consider the CC relaxation mechanism from Refs. |ll|, in which F in the 
action ( |2.1| ) is picked within the class of functions = F^{R,Q) of the form 

Here (3 and y are two parameters of the model, and n > and m > (usually taken to 
be integers) are two numbers characterizing a large class of functions that can realize the 
relaxation mechanism. To understand how the mechanism works in modified gravity, let 
us express the denominator of (|2.g|) in terms of H and q, 



B{H, q) = 2AH\q - )(g - 2) + [6 y H^l - qf . (2.10) 



Taking into account that each derivative of F in Eqs. ( |2.3D and (^^) increases the power 
of the denominator function B, it is easy to see that the induced energy density pF and 
pressure pF have at least one positive power of B in their denominators. Schematically, 
one can show that it takes on the form 

... flit) ^ f2{t) hit) 

with fi being time-dependent functions not containing the factor B. Notice that we need 
n < 2m in order that pF can increase with the expansion, i.e. when H is decreasing. In 
this case, the gravity modification F in the action ( |2.1| ) will be able to compensate the 
large initial vacuum energy density p\ in the field equations ( |2.2D , whatever it be its value 
and size, by pF = 2-E'o becoming large in a dynamical manner. It means that pF can take 
the necessary value to compensate for the big initial p\, this being achieved automatically 
because B becomes small during the radiation dominated regime ((/—)• 1), then also in 
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the matter dominated era {q — )■ and finally again in the current universe and in the 
asymptotic future, where H becomes very tiny. 

A few additional observations are in order. The fact that we search our relaxation func- 
tional among those of the form F{R,Q) is because they are better behaved. Indeed, let 
us recall that general modified gravity models of the form F{R, S,T), with R = g°'^Rab, 
S = RabR"^^ and T = Rahcd.R"'^'^'^ , are generally problematic as far as ghosts and other in- 
stabilities are concerned. This is because these theories introduce new degrees of freedom, 
which potentially lead to instabilities not present in general relativity (GR). They may e.g. 
suffer from the Ostrogradski instability, i.e. they may involve vacuum states of negative 
energy. In general they contain gravitational ghosts and can lead to various types of singu- 



larities. These issues are discussed e.g. in [21| and are reviewed in [22]. Fortunately, some 



of these problems can be avoided by specializing to F{R, Q) functionals involving only the 
Ricci scalar and the Gauf^-Bonnet invariant [^]. The reason is that all functionals of the 
form F{R,T — 4:S) are ghost free |24], a property which is shared by our functionals because 



g = R'^ - AS + T. 

Finally, the F{R,g) theories are thought to have a more reasonable behavior in the 
solar system limit [ p4| , p5[ . A detailed study of the astrophysical consequences of our model 
was presented in It shows that a huge cosmological constant can also be relaxed in 
astrophysical systems such as the solar system and the galactic domain. These studies 
were performed by solving the field equations in the Schwarzschild-de Sitter metric in the 
presence of an additional 1 /R term. The latter operates the relaxation of the huge CC at 
astrophysical scales in a similar way as the ~ 1/B term in equation ( p. 91 ) does it in the 
cosmological domain. One can show that none of these terms has any significant influence 
in the region of dominance of the other. Therefore a huge CC can be reduced both in 
the local astrophysical scales and in the cosmological one in the large. The solution in 
the Schwarzschild-de Sitter metric merges asymptotically with the cosmological solution, 
where we recover the original framework presented in Q. Furthermore, it was shown in 
jl2U that there are no additional (long range) macroscopic forces that could correct in 
a measurable way the standard Newton's law. This is a consequence of the dynamical 
relaxation mechanism and is in contrast to ordinary modified gravity models. The upshot 
is that the implementation of the relaxation mechanism should not be in conflict with local 
gravitational experiments. At the same time, we expect small deviations from GR only at 
scales of a few hundred kpc at least, which is in accordance with other authors |2^. As 
this scale is much larger than the scale where star systems can be tested at the level of 
GR, the physics of these local astrophysical objects should not be affected. 

2.1 Physical scales for the relaxation mechanism 

Remarkably, in order to insure that \p\ + Pf\ ^ \P\\ holds at all times, we need not fine 
tune the value of any parameter of the model. It is only necessary that the parameter /3 in 
the class of invariants ( |2.9D has the right order of magnitude and sign. This requirement 
has nothing to do with fine tuning, as fine tuning means to adjust by hand the value of the 
parameter f3 to some absurd precision such that both (big) terms p\ and pp almost cancel 
each other in (p.S]). This is not necessary at all for our mechanism to work because it is 
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dynamical, meaning that it is the evolution itself of the universe through the generalized 
Friedmann equation ( |2.8| ) what drives the F-term density pp towards almost canceling 
the huge initial p\ in each relevant stage of the cosmic evolution (i.e. in the radiation 
dominated, matter dominated and dark energy dominated epochs). In this way we can 
achieve the natural relation 



PAeff _ Pa + Pf{H* 



Pc 



Pi 



O 



p1 

p\ 



« 1, 



(2.12) 



for some value of sufficiently close to Hq. The smallness of the observed />Aefr as 
compared to the huge initial p\ thus follows from the right-hand side of ( 2.12| ) being 
suppressed by the ratio of the present critical density p\ = 3Hq / (SttG) as compared to its 
initial value in the early universe, which was of course of order p'^j^. By working out the 
explicit structure of (2.11) from (p.3[), and using the fact that H = —{q + 1)H'^, we can 
easily convince ourselves that all terms end up roughly in the form pp{H) ~ 13 H'^'^~'^^ . 
Although we are omitting here other terms, we adopt provisionally that expression for the 
sake of simplicity (see the next section for more details). Within this simplified setup, it 



follows that the value of H^: that solves equation ( 2.12 ) - which, as we said, should be close 
enough to the current value oi H ~ is approximately given by 



Hi ~ 



l/(2m-n) 



(2.13) 



We see that for n < 2m this mechanism provides also a natural explanation for the current 
value of H being so small: the reason simply being that p\ is very large! Moreover, we 
can make of order of the measured Hq provided (3 has the right order of magnitude 
(without operating any fine tuning). Notice that since the power mass dimension of /? for 



the models is 



and p\ 



|/3| = A^4-2n+4m ^ 

Mi, it follows that = Hq is related to Ai through 



M 

Mx 



2/(2m-n) 



(2.14) 



(2.15) 



For the simplest model this implies AA ~ \/MxHq, and therefore if Mx is near Mp 
we obtain Ai around the characteristic meV scale associated to the current CC density: 
10~^ eV. On the other hand, for the models F^ and we find in both 



niA 



,o^V4 



ipl) 



cases M ~ {Mj.HoY^^, which implies ~ 0.1 - 100 MeV for Mx in the ballpark of 
the GUT scale (~ GeV) up to the Planck mass Mp 
possible mass scale in the SM of Particle Physics. 



10 GeV. This is certainly a 



Notice that equations ( |2.15 ) and ( 2.15 ) have been derived under the assumption n ^ 
2m. Therefore, the conclusion that the current value of H is small because p\ is large 
cannot be inferred from equation ( p. 13 ) if n = 2m. Still, we shall see in the next next 
section that a more accurate treatment of the n = 2m models leads once more to the 
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conclusion that the current H is smah. What else can be learnt from equation ( 2.15| )? 
Notice that it can be rewritten as 

\n(Mx/M) 2(s-r) 

2m -n- 2 ^ \ / ^ 7^, (2.16) 

ln(7W/i7o) r + 42 ' ^ ' 

with Mx = 10" GeV, 7W = W GeV and i^o ~ 10""^^ GeV (for some integers r and s). 
Now, being Mx ^ Mp ^ 10^*^ GeV the GUT scale, we must have r ^ s ^ 19. For r ^ 0, 
it follows from (2.16) that ^ 2m — n < 1. But let us keep in mind that we could have 



r < 0. Since, however, the scale M. should not be too tiny^, it is natural to assume that it 
is not much smaller than the typical meV scale associated to the CC density, i.e. we must 
have M > mA ~ meV= 10~^^ GeV, which enforces r > —12. Thus e.g. in the extreme 
case s = 19 and r = — 12, we have 2m — n < 3. The upshot of this (order of magnitude) 
consideration is that, with r bounded in the approximate interval —12 < r ^ s, we should 
expect in general that the following combined "natural relaxation condition" is fulfilled: 

^ 2m-n < 0(1), (2.17) 

in which both naturalness of the mass scales and relaxation of the vacuum energy are 
insured. In other words, the above argument suggests that if n and m are taken as integers, 
we cannot assume large values for them (unless n = 2m) because the relation ( 2.17| ) could 



not be fulfilled. This is of course a welcome feature because it suggests that only the 
canonical cases are natural candidates, namely (n, m) = (0, 1), (1, 1), (2, 1), (3, 2) and (2, 2) 
at most - the last case of this series being the one which approaches the most to the 
upper bound imposed by the relation (2.17), but still the difference is of 0{1). The next 



model in the list, {n,m) = (2,3), tenses up a bit too much perhaps that bound, so for 
definiteness we stop the number of candidate models to the first five in the list. Of course 
the models with n = 2m for arbitrary n would do, but again it is natural to focus on only 
the canonical representative (n, m) = (2, 1) of this class. In general the class n = 2m is 



special because, as we can see from equation (|2.16D , it requires ^A = Mx- Therefore, if 
Mx ~ 10^^ GeV is a typical large GUT scale, then Ai must coincide with it. We may think 
of the n = 2m models as the class of "no-scale" relaxation models inasmuch as they do not 
introduce any other new scale beyond the one associated to the initial vacuum energy itself, 
Mx = (p\)^^^. Also interesting about the "no-scale models" is to note that if there is no 
GUT scale above the SM of Particle Physics, i.e. if Mx turns out to be just the electroweak 
scale Mw = O(IOO) GeV, then the relaxation mass parameter M will naturally take on 
the order of magnitude value Mw, characteristic of the electroweak gauge boson masses, 
for all the relaxation models n = 2m. Remarkably, it is the only situation where we can 
accommodate the electroweak scale Mw in the relaxation mechanism. 

To summarize, even though we have in principle five canonical candidate models sa- 



tisfying the natural relaxation condition (2.17), on the whole they involve only two new 



*If A4 could be very small, we would stumble upon the same tiny mass problem that afflicts e.g. all 
quintessence-like models Q], where mass scales as small as Ho ^ 10"'^'' eV are a common place. This is one, 
but certainly not the only one, of the most serious drawbacks of these models, another one being of course 
the need of extreme fine-tuning. In our case, however, we aim precisely at avoiding both of these severe 
problems. 
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physical scales for the dynamical adjustment of the cosmological constant, to wit: the two 
models ~ F2 are both linked to the same mass scale A4i close to the one inherent 
to the current CC density, A^i ~ m\ ~ 10~^ eV, whereas the two models ~ F2 are 
both associated to the scale A^2 ~ 0.1 — 100 MeV, which is a typical Particle Physics 
scale within the SM of electroweak and strong interactions. Thus both scales quite natural 
ones. On the other hand the fifth model Ff (and, for that matter, the entire F^ class) 
is a "no-scale" model which is able to operate the relaxation mechanism by using the very 
same scale as the one associated to p\ ~ M^, irrespective of Mx being a huge GUT scale 
or just the electroweak scale of the SM. 

As for the dimensional parameter y in ( 2.10| ), let us note that it is basically irrelevant 



for the present discussion. In Ref. Q, this parameter played a role only in the radiation 
dominated epoch and did not suffer any fine tuning either; it just had to take a value within 
order of magnitude. In practice the whole structure {y i?)^ of the last term in equation 



(2.10) was motivated by considerations related to getting a smooth transition from the 
radiation to the matter dominated epochs. Since, however, we are now comparing the 
class of models with the observations, all the relevant data to which we can have 

some access belong to the matter dominated epoch until the present time (apart from a 
correction from the radiation term in the case of the CMB data, which has been duly taken 
into account, see Sect. |3|). We have indeed confirmed numerically that the last term on 
the r.h.s. of equation (|2.10|) is unimportant for the present analysis. 

2.2 A closer look to the general structure of the induced DE density pp 

After we have presented the simplest version of the relaxation mechanism, it seems appro- 
priate to discuss a bit further the general behavior of the F-density pp order to better 
understand the relaxation mechanism in the different cases. Following the considerations 



exposed at the end of the last section, hereafter we will set y = in equation (2.10). We 



start by considering the the general structure of the -F-density pp = arbitrary 



FJ^ model, which as we know can be cast as in equation (2.11). However, by dimensional 



analysis and the explicit structure of ( |2.3| ) and ( 2.10| ), it is not difficult to convince oneself 



that it can eventually be brought into the more specific form: 



PF{H,q,q/H) 



(2.18) 



Here f{q) and g(q) are functions of q which are different for different -F^ models, but 
in all cases they are rational functions of q with a (multiple) pole at q = 1/2, specifically 
f{q) oc g{q) oc where for convenience we have defined the following expression 

that appears in the calculation: 

6:=2(g-i)(g-2). (2.19) 

For specific realizations of equation ( |2.18| ), see Appendix ^ Let us note that because of 
the pole at q = 1/2 all these models have an unstable fixed point in the matter dominated 
epoch, which is responsible for the universe to approach g — 1/2 for a long while during 
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that epoch and at the same time enhances dynamicahy the value of pp. As a resuh the huge 
value of the vacuum energy can be relaxed from p\ to the effective tiny |pAeff| ^ I Pa I' which 
can be identified with pAeff = Pa present time. When that pole is left behind during 

the cosmic evolution, the relaxation mechanism can still work effectively provided H tends 
to a very small value in the present universe. Notice that there are two terms in equation 
(2.18) that cooperate to fulfill this end: one of them is the overall ~ i^ff4ni-2n fg^^^^Qj.^ 
which contributes to the relaxation mechanism provided n is strictly smaller than 2m, as 



indeed required by the relation ( 2.17 ); and the other is the term q/H in the parentheses. 
For all the models such that n < 2m, the first term suffices and in this way we recover the 
result ( |2.13| ) which we had sketched in the simplified exposition of Sect. 2.1 . 

However, if n = 2m the simplified argument of Sect. p.l| cannot be applied. The 
class of these "no-scale" functionals is of the form ~ M^i?^™/i?™, where the relevant 
parameter is Mx, the same one as that associated to p\ ~ -^x- As previously noticed, 
this is the only case where the scale M defined in ( |2.14| ) could be identical to the initial 
vacuum energy scale Mx = {pa)^^^ ^ 10^^ GeV of the GUT at the early universe, which 
is a very interesting possibility in that here the mechanism for canceling the large vacuum 
energy density of the early universe would be naturally dealt with by the very same scale 
pertaining to the GUT transition at that time. For the n = 2m models, the factor q/H 



inside the parenthesis in equation ( 2.18 ) takes its turn in the relaxation mechanism. Noting 
that now equation (2.18) simplifies to 



PFiq,q/H)=f3 



{n = 2m) 



(2.20) 



it follows that the relaxation condition ( |2.12| ) for the present universe enforces H^, to take 
the value 



— m)q 

Pa 



\9{q)q\ ~ 1^1 ■ 
= 2m models, we have (3 



(2.21) 



M'x, 

and therefore the two factors cancel 



In this equation we have used the fact that, for the n 
which is of the same order of magnitude as p\ 
out approximately. Finally, g{q) is a dimensionless function, basically a number which 
can be taken of order 1 in this kind of consideration. Therefore, we conclude that the 
predicted value of the current expansion rate for these models is H^, ~ \q\, and therefore 
it can naturally be of order Hq. So, again, the relaxation mechanism works and requires a 
very small value for the present Hubble rate. 



3. Likelihood analysis from CMB, SNIa, BAO and H{z) data 

In the following we briefiy present some details of the statistical method and on the observa- 
tional samples and data statistical analysis that will be adopted to constrain the relaxation 
models presented in the previous section. First of all, we use the Constitution set of 397 
type la supernovae of Hicken et al. [27|. In order to avoid possible problems related with 
the local bulk flow, we use a subsample of 366 SNIa, excluding those with z < 0.02. The 
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corresponding XsNia function, to be minimized, is: 



XsNIa(p) 



366 

E 

1=1 



/^th(ai,p) - ^obs(ai 



(3.1) 



where = (1 + Zi)~^ is the observed scale factor of the Universe for each data point and 
Zi the corresponding (measured) redshift. The fitted quantity fj, is the distance modulus, 
defined asfi = m — M = 5 log + 25, in which dL{a, p) is the luminosity distance: 



di(a,p) 



da' 



c(l + z) 



dz' 



(3.2) 



with c the speed of light and p = {i^m, qo, qo/Ho) a vector containing the cosmological pa- 
rameters of our model that we wish to fit for. In equation (^.1|), the theoretically calculated 
distance modulus /^th for each point follows from using (|3.2D , in which the Hubble function 



is given by the generalized Friedmann's equation (2.8). Finally, /iobs 

(oj) and (Tj stand for 

the measured distance modulus and the corresponding la uncertainty for each SNIa data 
point, respectively. The previous formula ( |3.2[ ) for the luminosity distance applies only for 
spatially flat universes, which we are assuming throughout. 

On the other hand, a very interesting geometrical probe of dark energy is provided by 
the measures of H(z) [ po| l from the differential ages of passively evolving galaxies [hereafter 
H(z) data] ^. This sample of galaxies contains 11 entries spanning a redshift range of 
< z < 2, and the corresponding xii function can be written as: 



Xh 



11 

E 



Hthizi,p) - Hohsizi 



(3.3) 



Note that the latter set of data became recently interesting and competitive for constraining 
dark energy, see e.g. |28|. 

In addition to the SNIa and H{z) data, we also consider the BAO scale produced in 
the last scattering surface by the competition between the pressure of the coupled baryon- 
photon fluid and gravity. The resulting acoustic waves leave (in the course of the evolution) 
an overdensity signature at certain length scales of the matter distribution. Evidence of 
this excess has been found in the clustering properties of the SDSS galaxies (see |19], [^ ]) 
and it provides a "standard ruler" that we can employ to constrain dark energy models. 
In this work we use the measurement derived by Eisenstein et al. [19|. In particular, we 
utilize the following estimator 

2/3 

^(P) 



E{z) 



(3.4) 



[z2E(z,)]l/3 

with E{z) = H{z)/Hq the normalized Hubble rate. The previous estimator is measured 
from the SDSS data to be ^ = 0.469 ± 0.017, where Zs = 0.35 [or = (1 + Zs)'^ ~ 0.75]. 
Therefore, the corresponding Xbao function can be written as: 

[^(p) - 0.469]2 



Xbao(p) 



0.0172 

^For the current value of the Hubble constant, we use Ho = 100 ft = 73km/s/Mpc |28[, 



(3.5) 
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Finally, a very accurate and deep geometrical probe of dark energy is the angular scale 
of the sound horizon at the last scattering surface, as encoded in the location if^ of the first 
peak of the Cosmic Microwave Background (CMB) temperature perturbation spectrum. 
This probe is described by the CMB shift parameter |^, defined as: 



R=V^J ^-— = y^J -—. (3.6) 




The measured shift parameter according to the WMAP 7-years data [T| is i? = 1.726ib0.018 
at the redshift of the last scattering surface: zis = 1091.36 [or ais = (1 + zis)^^ ~ 9.154 x 
10~^]. In this case, the Xcmh function is given by: 

2 . . [Rjp) - 1-726]^ 
Xcmb(P) = • (3-7) 

Let us note that when dealing with the CMB shift parameter we have to include both 
the matter and radiation terms in the total normalized matter density entering the E(z) 



function in (3.6): 

nm{z) = nl{i + zf + n%{i + z)\ (3.8) 

Indeed, we have = (1 + 0.227A'^^) f^JJ, with N^, the number of neutrino species and 
Q^h'^ ~ 2.47 X 10~^. Therefore, at zis = 1091.36 and for three light neutrino species 
the radiation contribution amounts to some 24% of the total energy density associated to 
matter. For a detailed discussion of the shift parameter as a cosmological probe, see e.g. 

m. 

In order to place tighter constraints on the corresponding parameter space of our 
model, the probes described above must be combined through a joint likelihood analysis^, 
given by the product of the individual likelihoods according to: 

'^tot(p) = -^^SNIa X £h X -CbAO X £cmb ; (3.9) 

which translates in an addition for the joint function: 

X?ot(p) = XsNIa + Xh + XbAO + Xcmb • (3-10) 

In order to proceed with our minimization procedure, we would like to reduce as 
much as possible the free parameter space by imposing a prior, specifically we fix the value of 
the current mass parameter = Pm/ Pc- principle, is constrained by the maximum 
likelihood fit to the WMAP and SNIa data in the context of the concordance A cosmology, 
but in the spirit of the current work, we want to use measures which are completely 
independent of the dark energy component. An estimate of without conventional priors 
is not an easy task in observational cosmology. However, various authors, using mainly 
large scale structure studies, have attempted to put constraints to the parameter. In a 



^Likelihoods are normalized to their maximum values. In the present analysis we always report 2a 
uncertainties on the fitted parameters. Note also that the total number of data points used here is A'tot = 
379, while the associated degrees of freedom is: do f — Ntot — "-fit , where nflt is the model-dependent number 
of fitted parameters. 
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rather old paper Plionis et al. |33] using the motion of the Local Group with respect to the 
cosmic microwave background found ~ 0.30. Prom the analysis of the power spectrum, 
Sanchez et al. [34| obtain a value ~ 0.24. Moreover, [^] and [^] analyze the peculiar 
velocity field in the local Universe and obtain the values JlJn ~ 0->^0 ™d ~ 0.22 respectively. 
In addition, the authors of Ref. based on the cluster mass-to-light ratio, claim that 
lies in the interval 0.15 — 0.26 (see also |38] for a review). Therefore, there are strong 



independent indications for 0.2 < Q,^ < 0.3, and in order to compare our results with 
those of the flat ACDM we will restrict our present analysis to the choice = 0.27. If 
we fix the value of Q,^, then the corresponding vector p contains only two free parameters 
namely, {qo,qo/HQ). Note that we sample qo G [—1.24,0.1] and (Jo/Hq £ [—5,0.2] in steps 
of 0.001. 

4. Identifying the best relaxation models from observation 

In order to select the optimal relaxation models F^, defined in Eq. (p. 9]), from the phe- 
nomenological point of view, we are going to perform a likelihood analysis along the lines 
described in the previous section. Specifically, we will use a two-parameter fit of our mod- 
els. Namely, for any given model with fixed m and n, we have to look for the best fit 
values for qo and qo/HQ. 

4.1 Numerical solution of the models 

Remember that the parameters go and % are necessary in order to establish the initial 
conditions for solving the generalized Priedmann equation ( |2.8[) , in which pAefr is a compli- 
cated function of the form p\(,g{H, q, q/H) = p\ + Pf{H, q, q/H). Using H = —H'^{q + 1) 
one can see that p\eS = Pkes{H,H ,H) and hence the generalized Priedmann equation, 
despite its innocent appearance, becomes a third order differential equation in the scale 
factor a = a{t). Therefore, since the current is known, we need to input go and go for 
any given relaxation model F^. Notice also that the initial values of go and go must be 
consistent with the current value of which in our case is identified with 

OO _ pi + pF{Ho,qo,qo)/Ho _ . .... 

Pc 



for fiat space cosmology. Por the reasons indicated in Sect. 2.1, we will limit ourselves to 
analyze the five canonical models F^, F^, F2,Fi and F2 using the combined likelihood 
method described in the previous section. We present a sample of the likelihood contours 
in the {qo,qo/ Hq) plane for the individual sets of data on SNIa, CMB shift parameter, 
BAO and H{z) in Pigs. 1-3, whereas in Pig. 4 we display the combined likelihood contours 
for all the five models. A numerical summary of the statistical analysis for these models is 
shown in Table 1. In the next section, we provide more details of this analysis and discuss 
the obtained results. 



^Although equation (2.4) indicates that, in principle, the field equations are of fourth order, the self- 



conservation of pF, see (2.7) (and hence also of the effective vacuum energy pAcft), enables us to reduce the 



order of the field equations by one unit. 



- 13 - 



Table 1: Results of the overall likelihood function analysis, equation ( 3.10 ). The 1*** column 
indicates the various models. The ACDM model is included in the first row for comparison, 



although in this case we present the theoretical prediction based on equations (4.2) and (4.3). The 
2nd^ 3rd g^j^j ^th columns show the best fit parameters and the reduced Xtot- We use the prior 
fi^ — 0.27, fiat space cosmology and the (arbitrarily chosen) initial vacuum energy such that 
Pa/Pc — —10^'' for all the models. In the final column one can find various line types appearing in 
Fig. 4 where all models are plotted together. 



Model 


qo 


qo/Ho 


X?ot/377 


Symbols 


ACDM 


-0.595 


-0.887 






F? 


-0.506 ± 0.04 


-0.43 ± 0.03 


1.188 


red dashed 


Fl 


-0.500 ±0.03 


-0.36 ± 0.07 


1.187 


magenta long-dashed 


Fi 


-0.551 ±0.04 


-0.63 ± 0.07 


1.186 


black solid 


F? 


-0.678 ± 0.04 


-2.27 ±0.34 


1.188 


green dotted 


Fi 


-0.520 ±0.06 


-0.480 ± 0.05 


1.187 


black dotted area 



4.2 The Statistical Results 

In the upper left panel of Fig. |l| we present the results of our analysis for the model 
in the {qo^qo/ Hq) plane. The individual contours for each observable are plotted only 
for the la and 2a confidence levels in order to avoid confusion. In particular, the SNIa- 
based results indicated by thin solid lines, the H{z) results by thick dot-dashed lines, the 
BAO results by dotted lines and those based on the CMB shift parameter by thin dashed 
lines. The remaining panels show the statistical results for SNIa/BAO, SNIa/i!f(z) and 
SNIa/CMB. Using the SNIa/BAO data alone it is evident that the qo/Ho parameter is 
unconstrained within 2a errors. However, within la errors we can put some constraints 
(the best fit values are qo ~ —0.56 and qo/Ho ~ —0.64). On the other hand, utilizing the 
SNIa/i?(2;) data the best fit parameters are partially constrained within 2a: qo — —0.516 
and qo/Ho — —0.47. As can be seen in the lower right panel of Fig. |^, the above degeneracy 
is broken when using the SNIa/CMB data and practically the best fit parameters coincide 
with those of joint likelihood analysis, involving all the cosmological data. Indeed, for the 
F^ model we find that the overall likelihood function peaks at qo = —0.506 ± 0.04 and 
qo/Ho ~ -0.43 ± 0.03 with Xtot - 447.83 for 377 degrees of freedom. 

In the case of the "no-scale" Ff model we find that only the combined SNIa/CMB data 
can put constraints on the free parameters (see Fig. ^) , while the overall likelihood analysis 
provides qo = -0.678±0.04 and qo/Ho = -2.27±0.34 with xlot - 447.83 for 377 degrees of 
freedom. Concerning the Fl model (see Fig. |3|) we find that the corresponding statistical 
results are in a very good agreement to those of the F^ model. Indeed, the SNIa/BAO 
data put some constraints within la errors, go — —0.56 and qo/Ho — —0.66 while using 
the SNIa/-ff(2;) data we find qo — —0.530 and qo/Ho — —0.441. Again we observe that 
the joint likelihood function (mainly due to SNIa/CMB) peaks at qo = —0.50 ± 0.03 and 
qo/Ho —0.36 ± 0.07 with Xtot — 447.5 for 377 degrees of freedom. As for the F^ model, 
again we find that the comparison between SNIa/BAO does not place constraints on the free 
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-1 -0.5 -1 -0.5 

Figure 1: model: Likelihood contours in the {qo,qQ/Ho) plane. The contours are plotted where 
— 21n(£/£„iax) is equal to 2.32 and 6.16, corresponding to la and 2a confidence level. We do not 
plot the 3tT contour in order to avoid confusion. In the upper left panel we present the likelihood 
contours that correspond to the SNIa (solid lines), H{z) (yellow thick dot-dashed lines), CMB/shift 
parameter (green dashed lines) and BAOs observational data (red dotted lines). The remaining 
panels show the statistical results for different pairs. 



parameters while the SNla/ H{z) data put some constrains (even within Ifi): qq ~ —0.547 
and Qo/Hq ~ —0.63. As before, the free parameters of the model are tightly constrained 
by the SNIa/CMB data. The joint likelihood function peaks at qq = —0.551 it 0.04 and 
Qo/Hq = —0.63 ± 0.07 with Xtot — 447.12 for 377 degrees of freedom. To this end for the 
-F| model we find that the SNla/ H{z) comparison implies qo ~ —0.540 and Qo/Hq ~ —0.58 
while using the joint likelihood analysis we find qo = — 0.52ib0.06 and (Jq/Hq = — 0.48ib0.05 
with xiot — 447.53 for 377 degrees of freedom. Although we do not present individual 
likelihood contours for the and models, we can see their overall likelihood contours 
in Fig. ^ together with those of the other models. 

The summarized analysis of all the five canonical relaxation models is presented in 
Table |l] and in Fig. |^. In this combined figure we plot the la, 2a, and also the 3a, overall 
likelihood contours for the five the models under consideration. In it we can see a compact 
presentation of all our statistical results including their comparison with the corresponding 
values of {qoA^QoA/Ho) for the concordance or traditional ACDM model. Let us note that 
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Figure 2: model (or "no-scale" relaxation model): The Likelihood contours in the {qo,qo/Ho) 
plane. The different observational data are represented by different line types (see caption of Fig. 1 
for definitions). 



the (qAiz), qj^{z)/HQ) pair for the concordance model can be computed from the formulae: 

qK{z) = l^m{z) - 1 ; ^ = -I Ej,{z){l + z) = -\Et.{z)^^{z) [1 - ^^{z)\ , 

I Hq I az I 

(4.2) 

where 

^^™(-) = ^ = ; E^{z) = ^ = ^^l,{l + zf + l-^l,. (4.3) 

Pc{z) EX[z) Ho 

Therefore, for il^ = 0.27 we obtain (goAi Qoa/Hq) ~ (—0.595, —0.887), as quoted in the first 
row of Table 1. Let us mention that we have checked that using the earlier SNIa results 
(UNION) of Kowalski et al. does not change significantly the previously presented 
constraints. Finally, let us clarify that the chosen initial value p\ = — lO^'^ for the vacuum 
energy is completely arbitrary and the results do not depend on it because the relaxation 
mechanism is dynamical and hence works for any numerical choice of p\. However, in order 
to avoid instabilities in the lengthy numerical analysis involved in the solution of these 
models (in which the large quantity p\ almost cancels against the dynamically generated 
numerical value of pp) it is convenient not to choose p\ exceedingly large, but apart from 
this proviso any other arbitrarily selected value would do, as we have checked. 



o 
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SNIa/CMB ^ 



-1 



-0.5 



-0.5 



Figure 3: i^^ 'model: The Likelihood contours in the {qo,qo/ Hq) plane. The different observational 
data are represented by different line types (see caption of Fig. 1 for definitions) . 



5. The expansion history /future for the FJ^ models 

Since the current main cosmological quantities (scale factor, Hubble function etc.) exhibit 
a complicated scaling with the redshift, the absorption of the extra effects from the re- 
laxation model into an "effective dark energy" contribution, with a non-trivial EoS of the 
form uJes{z) = Pacs{z)/ p\efi{z) is possible with the effective pressure pacs{^) and energy 
density PAes{z) taken from Eqs. (|2.5D . The corresponding EoS for the matter component is 
not affected by the presence of the F-term, so that the behavior of the matter epoch is the 
expected one. Indeed, during the background evolution in e.g. the non-relativistic matter 
dominated era, the deceleration parameter q changes only very slightly to compensate the 
decreasing Hubble rate H ~ t^^, but q always stays around the value ^. Therefore, the uni- 
verse expands like a matter dominated cosmos despite \pW » pm- Analogously, q will vary 
minimally around the value 1 to ensure |/9Aeff| ^ I/^aI during the radiation regime. Eventu- 
ally, in the asymptotic future the relaxation of the CC is guaranteed by the smallness of H 
in the function B ( p. 10 ) at late times. The correct temporal sequence of the three cosmic 
epochs follows from the different powers of H in the function B, where higher powers are 
more relevant at earlier times: ~ for the radiation dominated epoch, and ~ for the 
matter and vacuum dominated epochs, although for the latter H becomes smaller than in 
the former because q is no longer close to 1/2. As we can see, all these dynamical features 
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Figure 4: The overall Likelihood contours in the {qQ^q^/ Hq) plane, computed from (3.9). The 
contours are plotted where — 21n(£/£niax) is equal to 2.32, 6.16 and 11.83, corresponding to Icr, 
2(7 and 3cr confidence level. The different cosmological models are represented by the following 
types: red dashed line, i^f black solid line, Ff green dotted line and F| black dotted area (see 
Table I). For comparison, the cross corresponds to the traditional ACDM model. 



are encoded in the structure of ( 2.10| ). The precise behavior of the relevant quantities in 



the various epochs becomes transparent in the various numerical examples considered in 
this section (cf. Figs. ^and|^, on which we will elaborate further below). 

In general, it is well known that one can express the effective dark energy EoS pa- 
rameter in terms of the Hubble rate, H{z). This function of the cosmological redshift z 
becomes known (numerically) after we explicitly solve the model as indicated in Sect.§. In 
the present case the structure of the effective EoS of the DE is quite cumbersome. First of 
all, from the formulae of Sect.|2| one can show that 

OJesiz) = 7-r = -1 + [1 + OJf{z)\ ^ = -1 + 2 , nrpQt \ ' (^"l) 



PKes{z) Pi + PF{Z) Pa + '^K 



z 



where 



In view of equations ( |2.3D and (|2.4| ) , the previous expressions are complicated functions 
of z. From { f).l\ ) and ( ^.2| ) it is patent that if we would have E'^q = (1/3) E^-, then the EoS 



,,PFiz)_ 1 E\iz) 
Pf{z) 3 E'>q{z) 
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Figure 5: The effective equation of state UcS as a function of tlie redshift z for the models 
with initial conditions taken from Table 1 and with Q.'^-^ = 0.27. The horizontal line lOcB = ^1 
corresponds to the ACDM model. 



of the induced DE density would he cop = —1, and at the same time the EoS of the 
measurable effective vacuum energy would also be cjeff = —1- However, as we can see 
from equations ( |2.3| ) and (p.4[), the relation E^q = (1/3) £'*j is only satisfied by the F-terms 
on the respective r.h.s. of these equations, but not by the remaining terms. Therefore, 
the effective EoS of the measurable DE density is expected to have a non-trivial behavior, 
and we must check if this behavior is still sufficiently close to the ACDM result u}\ = — 1, 
particularly at z = 0, and also determine what is its behavior as a function of the redshift. 

In order to visualize the redshift dependence of the effective EoS parameter, we compare 
m Fig. I the various cosmological models indicated in Table I of Sect.|^. One can divide 
the evolution of the cosmic expansion history in different phases on the basis of the varying 
behavior of the and ACDM models. We will investigate such variations in terms of 
the deceleration parameter, q{z) = —1 + d{ln H) / d{ln (1 + z)), which is plotted in Fig. |6|. 
Below we briefly present the cosmic expansion history of the models studied here. Note 
that we always compare with the concordance ACDM cosmology. 

Although we do not focus here on the details of the radiation dominated period within 
our model, let us recall that it has been discussed in Sect. 5. 2 of Ref. [|ll|. From the ap- 
proximate effective equation of state of the radiation epoch - see eq.(5.8) of that reference 
- one can see that it behaves radiation-like, i.e. Wefr — 1/3, and smoothly connects with 
(Jeff — in the matter dominated epoch during equality. Therefore the behavior is perfectly 
compatible with that of the ACDM model. This is corroborated by the numerical solution 
of the full field equations [^] . Finally, we point out that there is no significant effect that 
could alter the BBN phase of the radiation epoch because, as indicated previously, the 
form (p.g|)-( |2ri0 ) insures that the dynamical relaxation mechanism automatically reduces 
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Figure 6: The deceleration factor q(z) as a function of the redshift z for the ACDM model 
with = 0.27 and for the i^," models with initial conditions taken from Table 1. Note that 
for z > 0.2 the curves of the models ACDM and F2 are overlapping. 



the big initial CC to a tiny value at all epochs, starting from the radiation epoch, going 
through the matter epoch until our present epoch. The density of vacuum energy during 
BBN therefore was completely subdominant and could not have any measurable effect on 
the standard light element abundances. 

Prior to the radiation epoch we have the inflationary epoch. The latter should not 
be affected by the relaxation mechanism. This can be seen once more from the general 
structure of the functional ( p.9D -( p^ ). In the inflationary epoch, H is very large and q is 



near —1. Therefore the behavior of F in equation ( p.9[ ) is of the form F ~ l/H^"^ which 
after applying the natural relaxation condition (|2.17 ) implies it is severely suppressed as 
|F| < l/H^"", at least. Thus, the universe is then completely controlled by the physics 
of the inflationary epoch, whatever it be, without receiving any interference from our 
relaxation mechanism. This mechanism only starts working when the universe leaves the 
inflationary epoch and enters the radiation dominated one, since then q — t- +1 and this 
triggers the first large contribution from F that compensates for the value of the huge CC 
left over at the end of the inflationary period. Although we do not attempt to describe 
here how the inflationary epoch transited into the radiation dominated one, in the original 
formulation of the relaxation model [4] it was suggested that the functional F in equation 
( |2.9| ) could contain an additive polynomial A{R) of the Ricci scalar. In its simplest non- 
trivial form it would just entail a term proportional to R^, which would not alter the 
relaxation mechanism in the radiation epoch. The advantage of this addition is that it could 



connect this relaxation model with a Starobinsky's type of mechanism for inflation [4C] and 
subsequent modifications thereof, see e.g. [41| and references therein. 
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5.1 Cosmic acceleration and effective EoS analysis 

In the following, we remark some features that can be observed from the behavior of the 
effective EoS and deceleration parameter of the various relaxation models under consider- 
ation in Figs. |5| and ^, respectively: 

• model: For z > 1 the deceleration parameters (both for and ACDM) are 
positive with q(z) > Qa{z), which means that the cosmic expansion in the F^ model 
is more rapidly "decelerating" than in the ACDM case. At z ~ 0.8 both models are 
starting to accelerate. Between 0.2 ^ z < 0.8 the deceleration parameters are both 
negative with q{z) < qA{z). Now we will focus on the evolution of the effective EoS 
parameter. In particular, we find that at early enough times z > 3 the effective EoS 
parameter WeffC-z) of the F^ model approaches zero while we always have u}a{z) = — 1 
for the ACDM model. At z ~ 1.8 the F^ effective EoS parameter crosses the phantom 
divide line {ujes = — 1) and it stays there for some time (0.2 < z < 1.8). Close to the 
present epoch, the effective EoS parameter crosses the phantom divide again and it 
behaves quintessence-like at present, where tdesi^) ~ —0.92. In the future, it sustains 
this quintessence-like behavior and tends to uj^ ~ —0.83 . 

• Fi model: this case behaves qualitatively very similarly to the Fi model with a 
transient uj^s < —1 phase in the past and quintessence like behavior in the asymptotic 
future. It also behaves quintessence-like at z = 0, where it takes essentially the same 
value of a;eff(0) ^ —0.91 as in the previous model. The asymptotic EoS behavior is 

• ^2 model: for < z < 3 the deceleration parameter practically coincides with that 
of the concordance A cosmology, which means that the cosmic expansion in the F| 
model mimicks that of the ACDM case. Also, the effective EoS is not far from — 1 
in this region, but it stays above the phantom divide line. Specifically, the current 
value is u;eff(0) — —0.96. For z < (future) the deceleration parameters are both 
negative with q{z) > qA{z), which means that the cosmic expansion in the ACDM 
model is more rapidly "accelerating" than in the -F| case. The latter approaches a 
quintessence like future with ~ —0.73. 

• F^ model (or "no-scale" model): for 1 < z < 3 the deceleration parameter remains 
close to that of the concordance A cosmology. Then at z ~ 0.7 both models are 
starting to accelerate and between < z < 0.7 the deceleration parameters are both 
negative with q{z) > qA{z). Near the present epoch, the effective EoS parameter 
"^eff ~ — Ij in particular the current value is a;eff(0) ~ —1.08. Interestingly, in the 
future the F^ leads to an apparent singularity {iOes{z) <^ —1) caused by the pole 
that Weff has at z ~ 5.7 due to pAcS changing its sign. There is, however, no physical 
divergence related to it because all energy densities remain finite. The dramatic 
change of qualitative behavior of this model with respect to the others is nevertheless 

*This is no surprise, it is actually a common feature expected for general AXCDM models (as the 
relaxation models indeed are), see for details. 
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remarkable. In particular, the fact that behaves mildly phantom-like at present 
could explain the persistent tilt of the EoS data pointing slightly below — 1 at z = 0, 
if eventually confirmed by the observations. 

• F| model: close to the present time and in the future this model performs very 
similarly to the model with a;eff(0) — —0.93. However, the transient Wgff < — 1 
phase between z w 0.3 and z w 2.5 started earlier in the past, roughly at the same 
redshift where the model switched from the phantom- to the quintessence-like 
EoS. 

5.2 Phase space analysis for the "no-scale" models 

The "no-scale" models n = 2m, like the F^ discussed in the previous section, deserve a 
closer attention. From the peculiar structure of ( 2.20| ), we see that pp is a function of 



only two arguments: pp = pf{QtQ/H). This fact enables us to analyze the F^ models 
in a {q,q/H) phase space diagram. This is generally not the case for the other models, 
for which pp is in general a function of three (not just two) independent arguments pp = 
pp{H,q,q/H), see equation ( p. 18 ) and also Appendix If we focus once more on the 



canonical case F^, we explicitly obtain 

(5.3) 



pp{q,q/H) = ^ 



(g^ - V + -5q-l)- [Iq-" - 6g + 5)-i- 

Hb 



with /3 = M|- and b given by ( 2.19[ ). If we compare with (2.18), here we have f{q) 



6(g^ — + 9g^ — bq — l)/6^ and g{q) = — 6(2g'^ — 6(7 + 5)/fe^, and from the value of go taken 
by this model (see Table 1) we can see that f{qo) and g{q{)) are of order 1 for that matter. 
The corresponding phase space diagram is depicted in Fig. ^. The background evolution 
follows mainly from pp -\- p^^ = 0, which yields q/H as a function of only q, see Fig. |^ 

When q/H is very small, the deceleration does not change much. This happens 
around q ^ ^ corresponding to the matter era. For smaller values of (7 < ^, the time 
derivative q is negative and q will decrease with time. Finally, q will vanish again in an- 
other fixed point, which can be quintessence-like (g > —1) or phantom-like {q < —1), or 
q/H ^ —00 diverges at around q « —2.05. The type of final state depends on the value 
of p\/{2(3) (the values are attached to the curves in Fig. |^ or equivalently on the current 
values of qo and %. The smallness of the late-time Hubble rate Hq results from the stan- 
dard evolution of H ^ 1^^^ during the matter epoch, which corresponds to the unstable 
fixed point {q = ^, q = 0)- 



6. Conclusions 

In this paper, we have confronted the class of modified gravity models of the form F^(R, Q) 
- see equation (|2.9| ) - with the main observational data. We have shown that the set of 
all models of this kind with n < 2m are possible candidates capable to implement the 
dynamical relaxation of a large cosmological term, or vacuum energy density, irrespective 
of its value and size - typically p\ ~ (with Mx > 10-*^^ GeV). We have adopted the 
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Figure 7: Phase space diagram for tlie case n = 2, to = 1. We plot qjH as a function of 
the deceleration q for different values of p\/{2l3) (attached to the curves). The matter cosmos 
with q = \ (I'ight vertical dashed line) is an attractor in the past. The red blob denotes the 
best-fit point (go = —0.68, go = —2.3) found in the statistical analysis of Sect. 3. Moreover, we 
show q/H (dashed orange curve) for the ACDM model without radiation, which has also a stable 
fixed point at g = — 1 corresponding to de Sitter space-time (left vertical dashed line). In this plot 
q is decreasing/increasing with time in the lower/upper part divided by the horizontal dashed line. 



point of view that such a huge vacuum energy is an integral part of the energy budget 
of the universe, which was unavoidably deposited in it by the quantum fluctuations and 
the important phase transitions of the early times. All the usual modified gravity models 
we know of, and in fact most DE models of all kinds proposed in the literature, implicitly 
assume that this gigantic energy is just canceled by some extremely fine tuned counterterm 
such that the measured value is the tiny number ~ 10"'*'' GeV^ (in particle physics 
units), otherwise the cosmos could not evolve in the observed standard manner. Usually 
the measured value is then linked to some late-time effective gravity modification or to the 
residual value of some scalar field potential. However, the preposterous cancelation that 
was implicitly assumed to get rid of the huge initial vacuum energy density is generally 
recognized by theoretical physicists as completely unnatural and unacceptable. 

Quite in contrast, the relaxation mechanism that we have put to the observational test 
in this paper is able to counterbalance the large vacuum energy (whatever its value and size) 
in a totally dynamical way (hence without any sort of fine tuning) thanks to an effective 
action which is modified by the presence of the F^[R,Q) terms. By carefully comparing 
these relaxation models with the most recent observational data on SNIa, CMB, BAO and 
high redshift H{z), we have shown in particular that the simplest candidates in the list, 
namely F^, F/, F|, and F|, can provide a background cosmic evolution very close to 
the concordance ACDM model. In other words, these models yield a fairly standard-like 
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cosmological evolution going through the normal radiation dominated, matter dominated, 
and CC dominated epochs, they also exhibit an effective EoS behavior very near to — 1 
around z = 0, and finally (and this is of course the main point to be stressed here) they do 
all this without any need of enforcing fine tuning in order to get rid of the huge vacuum 
energy injected in the cosmos during the early stages of its evolution. In particular, the 
"no-scale" relaxation models, i.e. the entire class F^, can provide a natural relaxation of 
the vacuum energy without introducing any other energy scale except that of the vacuum 
energy itself. The other models, with n 7^ 2m, involve a mass scale which, for the 
canonical candidates mentioned above, is of the order of a typical Particle Physics scale in 
the SM of electroweak and strong interactions. Furthermore, in contrast to quintessence- 
like models, these relaxation models do not involve extremely tiny mass scales. Actually, 
the typical scales that are required are fixed in order of magnitude only, and therefore 
there is no need at all of exerting fine-tuning in any of the F^ models so as to dynamically 
reduce the huge initial value of the vacuum energy density p\ ~ into the very small one, 
~ 10"^'' GeV^, measured in our current universe. The class of models studied here can 
be thought of as a prototype for eventually solving the "old CC problem" without departing 
significantly from the expansion history of the concordance model. However, more realistic 
models are of course needed in order to provide a true solution to the cosmological constant 
problem from the point of view of fundamental physics. 
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A. Computing pp for the canonical models: , , , F2 , F2 

Here we quote the explicit expressions for pp = 2-E''o for our favorite F^ models in Table 



1 of Sect. 4.2, which are obtained from explicit computation of (p. 3D. To avoid too lengthy 



expressions we set y = in equation (2.1C) as it is not important for the matter dominated 



epoch or any time after it, and in this way they all take the general form pF{H,q,q/H) 
given in ( p. 181 ), with the notation ( p. 19 ). The explicit results read as follows: 

1. Case n = 0, m = 1, with = M^: 



(5(?2-3g-5)-(V_l0g + 7)-|^ 



2. Case n = 1, m = I, with \/3\ = 



/3 

PF 



4(4g^ - + 2q + A) + (V - 12?^ + 18g - 11)^7 
2 Ho 



(A.l) 



(A.2) 
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3. Case n = 2, m = 1, with |/3| = Ai^ (see Sect. 2^ for a detailed discussion of this 
special case): 



PF 



6/3 



{q^ - iq^ + - 5^ - 1) - (2^^ -6^ + 5) 



Hb 



4. Case n = 3, m = 2, with \/3\ = M^: 
3(<7-l)/3 



PF 



3{q - l)(V-6g2 + 5g + 6) 

+ (V - 12g^ + 28g2 - 36q + 19)-^ 



Hb 



5. Case n = 2, m = 2, with \/3\ = M^: 
/3 



PF 



8 63//4 



(309"^ - nq^ + 44g2 + 37g _ 38) 
- 2(12/ -48g^ + 



82g + 31) 



(A.3) 



(A.4) 



(A.5) 



In all cases is a mass scale associated to the relaxation mechanism, which depends on 



the particular model, see equation ( 2.14 ). The numerical value of that scale (fixed in order 
of magnitude only) is determined by the value of p\. The models listed above have been 



chosen because they all satisfy the "natural relaxation condition" ( 2.17 ), and for all of 
them Ai takes a characteristic value which ranges from a mass scale of the SM of Particle 
Physics to the GUT scale associated the big initial value of the vacuum energy, p\ 



See sections 2.1 and 2.2 for details. 
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